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00 ■ Abstract 



^ . In connection with numerous X-ray and neutron investigations of some high 

P^ ■ temperature superconductors (YBa2Cu306+x and related compounds) a non- 

^^ ■ trivial part of the structure factor, coming from partly disordered Cu-O-. . .- 

O I 0-Cu chain fragments, situated within basal planes, CuOx>, can be a subject of 

Qv^ I theoretical interest. Closely connected to such a diffusive part of the structure 



factor are the correlation lengths, which are also available in neutron and X- 

H ■ ray diffraction studies and depend on a degree of oxygen disorder in a basal 

i-^ ■ plane. The quantitative measure of such a disorder can be associated with 

t^ ' temperature of a sample anneal, Tg, at which oxygen in a basal plane remains 

J-^ ' frozen-in high temperature equilibrium after a fast quench of a sample to room 

or lower temperature. The structure factor evolution with x is vizualized in 

rS I figures after the numerical calculations. The theoretical approach employed 

^ in the paper has been developed for the orthorhombic state of YBCO. 
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I. INTRODUCTION. 

While the regular Cu02 planes are commonly believed to be responsible for both fasci- 
nating phenomena, the high-T^ superconductivity (HTS) and high-T/v antiferromagnetism, 
the oxygen deficient planes of the YBa2Cu306+a; compounds play an important role too, 
taking a control over the mentioned phenomena. HTS is associated with the orthorhombic 
state of YBCO which is consistent with existence of linear chain fragments, Cu-0-. . .-0-Cu, 
in oxygen deficient planes, CuO^;. Even in the case of twinned crystals the substance consists 
of macroscopical domains where chain fragments (CFs) are arranged in chains (conventional 
notation is 0(1) and 0(5) for oxygen positions in and across the chains, respectively). There 
are three orthorhombic phases, Ortho-I, II and III, which have been reliably identified ex- 
perimentally. Concerning the Ortho-II and Ortho-Ill phases, their identification in neutron 
and X-ray diffraction studies JIj-Q has been reported. The part of the structural phase 
diagram in the x-T plane is shown qualitatively in Fig.l (cf Refs ^-0). The T coordinate 
plays the role of temperature of oxygen equilibrium in the basal plane. It could be of the 
order of room temperature, Tr, if a cooling procedure were done slowly. Otherwise, it would 
be the so-called quenching temperature, Tq, if a sample were annealed for a long time and 
then quenched rapidly from Tg to room or lower temperature. The role of Tr is impor- 
tant: Oxygen atoms in deficient planes do not diffuse below Tr (cf, however, the elegant 
investigation by Jorgensen et al [Q). T^'s have varied up to 350°C in experiments by Veal 
and co-workers ( [^[101 and references therein). In order to prepare the Ortho-Ill phase the 
samples have been quenched at temperatures above 500°C 0. 

As mentioned above the constituents which form the oxygen deficient plane are the 
CFs. In the tetragonal phase they are rather short and could be assumed to be situated at 
random were it not for evidence of a superstructure at some low oxygen content, x ~ 0.35, as 
reported in |[TT| , p!2[] (cf [|13|, however). Most likely, the origin of the herring-bone arrangement 



of the very short Cu-O-Cu chain pieces is simply owing to the Coulomb interaction of such 



fragments which are simultaneously quadrupoles (see |I^). 

In the Ortho-I phase the distribution of CFs over their lengths is similar in all the chains. 
The single-cell structure is typical for the right part of the phase diagram (0.75< x <1 at 
Tq ~ Tr). Between Tetra and Ortho-I phases in Fig.l there is the area of the Ortho- 
II phase, which simply means a double-cell structure with oxygen-poor chains alternating 
oxygen-rich ones. However, the phase transition lines, restricting this area, are not strictly 
determined up to now: There is still a question about the type of the Tetra-to-Ortho-II 
phase transition. Either the system could undergo the phase separation, ie, the first order 
phase transition with phase co-existence in the form of the "gas" (Tetra)-" liquid" (Ortho-II) 
mixture or the phase transition through a growth of Ortho-II domains. More acknowledged 
is the first scenario, see, for instance, |P,|6|, p!5| , [l6| , but this conclusion has been done in a 



frame of a lattice-gas model, the so-called ASYNNNI model, which ignores intrinsic, ie, 
fermionic degrees of freedom of CFs. As for the second scenario the domains of four various 
kinds would be locally arranged in the double-cell structures, oriented along 6-/a-axis and 
shifted across the chains by a single period. Another boundary, confining the Ortho-II phase, 
would divide the regions of the Ortho-I and Ortho-II phases were it not for the possibility 
for the Ortho-Ill phase to exist between them. The periodicity of the latter corresponds 
to one oxygen-poor chain following two oxygen-rich chains. However, in order to search 
the Ortho-Ill phase experimentally the way of a sample preparation should be indirect. 
Actually, a slow cooling is accompanied with a formation of the Ortho-II phase, which locks 
the periodicity, making the kinetics of the phase transformation into the Ortho-Ill phase 
to be embarrassing. To avoid such a difficulty a slow cooling in the atmosphere of higher 
oxygen pressure at the first stage can be switched over to a fast quench from rather high 
temperature (>500°C) at the second stage Pj. There was a theoretical conjecture concerning 
a possibility for long wave-length modulated structures to exist in oxygen deficient planes 
(see [^ and references therein). However, the afore-said impossibility of oxygen migration 
at low enough temperatures {T<Tr) makes such hypothetical modulated structures the 
subject of an abstract interest. 

Note that the Ortho-I phase can be interpreted as a disordered Ortho-II phase with 
a short-range order. Certainly, the order parameter of the Ortho-I phase is of the higher 
symmetry than that of the Ortho-II phase. 

To reveal the regions of stability of Tetra, Ortho-I and Ortho-II phases and to construct a 
structural phase diagram theoretically a two-dimensional lattice-gas model, the ASYNNNI 
model, has been used as a prerequisite by de Fontaine and colleagues. The model pre- 
sented in Fig. 2 (on the parameters, entering the model, see, for instance, Pfljl^) invokes 
the nearest neighbor oxygen-oxygen repulsion (Vi, oxygen in the "corner" configuration), 
the next-nearest neighbor attraction (— V2, in-line configuration) and also the next-nearest 
repulsion without any Cu ion in between the corresponding oxygen atoms (V3, "vis-a-vis" 
configuration). As mentioned above the lattice-gas model ignores a fermionic origin of micro- 
scopical properties of CFs. Also this model is not the best candidate for a proper description 
of a structure formation: The superstructure at small x is beyond the ASYNNNI model 
predictions. Nevertheless, the lattice-gas approach is reasonable for the orthorhombic side of 
the phase diagram and it can be phenomenologically generalized to incorporate the fermionic 
degrees of freedom into the effective interaction of oxygen atoms within CFs. Generaliza- 
tion is based on the result obtained in the soft X-ray absorption spectroscopy investigation 



of untwinned YBa2Cu307 and reported in |I9|: The concentration of oxygen holes in the 



Cu-0 chains has been estimated around 30%. This experimental estimate is essential for 
applying the charge transfer mechanism discussed theoretically in ||2^ and partly based on 



the experimental investigations 0,0] of the monovalent Cu amount in YBa2Cu306+a;- A 



short resume of the mechanism proposed in [Q is as follows: Existence of the CF, involving 



n+1 Cu^+ and n O atoms, m of which are 0^~ and n—m atoms are 0~, is consistent with 
m—1 holes transferred from the fragment. In the limit of very long CFs {n^ 1), ie, the 
case investigated in Ref |jl9[|, the free energy per one oxygen site of the fragment, f{n,m), 



depends practically of m/n and has a minimum around ^ ~ 0.70. It allows to represent 
the free energy in the following form: 

/Km) = /„ + ^(^-en)' (1) 

where /„, a„ and ^„ go to the fixed values when n ^ oo. One employs Eq.(|l]) to get a 
reduced information about CFs, expressed in terms of their lengths n: 



e 



-n4>{n)/T ^ y^ g-n/(n,m)/T ^2) 

This procedure is convenient to build the statistical mechanics of the CF ensemble because 
the configurational entropy is associated with the lengths of CFs. The reduced free energy, 
0(n), reads at n^^l: 

T 27mT 4>2 T . . 

d){n) = fn m ~ 01 H In n (3) 

^^ ^ ■' 2n an n 2n ^ ' 

Expression (^ could be compared with the analogous one easily obtained from the lattice- 
gas model mentioned above. Actually, the attraction energy — V2(n — 1) plays the same role 
as n(j){n) in formation of the CF of length n. So, the analogue of Eq.(^) is 

4>{n) = -1^2 + - (4) 

n 

Deriving Eq.(|D we also include the energy loss, Vs-cu, due to the 3-fold coordinated Cu ions 
at the ends of a CF, which yields the evident change ¥2/11 -^ V^jn = (V2 -|- 2V3_cu)/'^- The 
first terms in the r.h.s. of Eqs.(^ and (^) are not relevant if oxygen atoms predominantly 
occupy 0(1) positions: For the orthorhombic state of YBCO, prepared at moderate (not 
too high!) Tq, 0(5) positions were proved |^ to be occupied poorly. This is the case for 
which one assumes the total amount of oxygen to be fixed, hence, the first, n-independent, 
term of Eqs.(|^)-(0) does not distinguish between various configurations of CFs. The second 
terms in Eqs.(^-(^ are important and their meaning is clear: CFs would show a tendency 
to increase their lengths were it not for the third term in Eq.(^) and, in addition, for the 
configurational entropy of the CF ensemble. Aforementioned temperature, T = Tq, is above, 
at least, than Tr. 

In order to complete the scheme presented by Eq. (|^) or (^) the following parts of the free 
energy must be involved into a theoretical consideration: First, the configurational energy 
term due to various configurations of CFs, and second, caused by the interaction of oxygen 



atoms in different chains, ie, the Va-term. Partly, this hne of attacking the problem has been 
realized in p4| , p5| and subjected to the experimental check by Boucherle et al [^]. Such a 



general theoretical approach as well as some important applications to the experimentally 
observable values, reported already and foreseen, will be presented in the following sections. 
There is also a more straightforward approach to the problem of CFs in the orthorhombic 
phase of YBCO. Instead of employing phenomenological schemes of Eq.(^ or (^) it would be 
possible to start with some microscopical model for an alternating Cu-0 CF, which should 
belong to a family of strongly correlated fermionic systems despite of its finite length. A 
good candidate could be the Emery model [^, or its Kondo-like version [Q, or the Id 



model of the Zhang- Rice type [^, derived from the original Emery model: All are well 



adapted to a proper description of the alternating CFs. However, to adjust such models to 
the reality one needs to invoke some additional short-range interactions, ie, the nearest Cu-0 
and 0-0 interactions, the bare energy of a charge transfer from chains into Cu02 planes 
and the energy loss because of the 3-fold coordinated Cu atoms at both ends of CFs. The 
problem seems to be hopeless because of too many energy parameters invoked in, however, 
the choice of the model parameters is based only on few relevant properties of CFs, the 
most important among them is a true position of the free energy minimum around m/n ~ 
0.7. The attempt to apply the Kondo-reduced version of the Emery model has been done 
in Ref |2^: In addition to the above-mentioned criterium the plateau- like behavior of nh vs 



x, imitating the 60 K Tc vs a; plateau in YBCO, has been employed as a crucial point (nh 
the concentration of holes leaving the chains). 

Certainly, a big advantage of the microscopical model approach is in a possibility of direct 
calculations of such physical quantities as rih, n-^mi the concentration of paramagnetically 
active CFs, studied by polarised neutrons |^ , magnetic correlations within CFs and so on. 



However, the chain, consisting of more than 9 copper and 8 oxygen atoms, are not available 
to be treated by means of the exact numerical diagonalization algorithm by Lanczos (see 
the review paper [^ for the technical details of the method). Nevertheless, the quantum 
Monte Carlo algorithm seems to be promissing in studying of microscopical properties of 
CFs at non-zero temperature. 

In this paper we mainly follow the line of the generalized lattice-gas model, whose pecu- 
liarities are expressed by Eq.(^. 

The principal results of the paper are collected in Section HI, while all the important 
spade-work is done in Section II and Appendixes. Section II starts with application of 
the linear programming method for selecting CFs of optimal lengths. But the method 
works reliably only at rather low temperatures when the contribution of the configurational 
entropy can be neglected. Then a more realistic situation is considered: The average lengths 
of chain and vacancy fragments are calculated as a function of oxygen filling in a chain vs 
temperature; the in-chain structure factor as well as the in-chain oxygen density- density 



correlation functions are analysed. In Appendix B we represent the equations of statistical 
mechanics of the ensemble of CFs arranged orthorhombically by means of the distribution 
of CFs over their lengths. In Section III the interchain interaction is incorporated by means 
of the high-temperature expansion of the free energy. The final results are partly expressed 
in the analytical form (Section III), partly visualized after numerical calculations. They are 
presented in Section IV, which also includes the set of plots for some observable physical 
quantities, playing the role of the reference set. 

Before starting the main part of the paper, let us introduce some definitions: N^, the 
total number of chain fragments, consisting of m+1 Cu atoms and m O atoms. It will be 
convenient, especially for further application to the statistical mechanics of the CF ensemble, 
to incorporate one of the vacant of oxygen Cu-Cu links, surrounding CF, to be certain, to 
the right, into it. These numbers and corresponding "probabilities" Xm = Nm/Af£ satisfy the 
evident constraints: 

1 

and 



— J2{m + l)Nm = E ('^ + 1)^- = 1 (5) 

J\lf 



S=—rJ2^rn=J2^rn = '^-X (6) 

where A/^ is the total number of Cu atoms in a chain. S is, evidently, the concentration of 
oxygen vacancies in a chain. Many experiments in YBCO have been performed at the fixed 
oxygen content: Because of a possibility for oxygen to fill the chains inhomogeneously one 
redefines Eq.(^ by taking a sum over the chains: 

^ = T^EE4"^ = i-^ (7) 

where a labels the chains and A/'ch is their total number. 



II. IN-CHAIN DISTRIBUTION OF CFs OVER THE LENGTHS. 

A. A hypothetical low-temperature case. 

Using the intrinsic free energy of a CF, Eq.(|D, let us elucidate an auxilary problem, 
namely, what kind of CFs would be preferable if the configurational entropy as well as the 
interchain interaction were neglected. The free energy of the ensemble of non-interacting 
CFs (per one copper site) can be expressed as a linear function of the set {xm}, which is 
also subjected to the /mear constraints (|§)-(|^): 

/ = E m(l){m)xm = E (" - rlnm)x„, (8) 

m=l m=l 



where m = 02, t = T/2. / can be minimized by making use the hnear programming 



method [0. In the physics of the CF ensemble this procedure would reflect the properties 
of non-interacting constituents, ie, CFs, if they were supposed to be annealed at rather low 
temperature. Note, however, that such a case has no direct bearing on the YBCO behavior, 
because oxygen of deficient planes would be supposed to be frozen below Tr. 

The linear programming theory says |^ that a function, / in this case, will attain its 



minimum on a unique vertex of the convex polytope V defined by constraints (§)-(|^). Any 
vertex [i,j] defined by Eqs.(|^)-(|^) can be represented by a pair of non-zero probabilities 
i^<3)■ 

(j + l)(l-x)-l l-(i + l)(l-x) ^ .^ .. .. ,^, 

Xi = ^ ^^ -^ , x, = . ^. -, x„ = if m ^ «, m ^ J (9) 

Certainly, not all the pairs are the candidates to be vertices, a true vertex of this particular 
problem is composed of two non-negative probabilities. This results in the following available 

set: 

1 1 — S" 

j> K>i, K = 1 = —— (10) 

1 — X D 

True vertex [i, j] defines /{«, j}, the function to be optimized, as follows: 

/{2,j} = (l-x)U-rf^lnj + ^lnzjj (11) 

It is proved in Appendix A that /{«, j} attains the global minimum on the manifold of 
vertices, defined by Eq.(D with i and j both non-zero, on unique vertex [/T;_,fi;_|_] {k± are 
the two nearest integers to k): 

/{k_, k+} = {\ — x)\u — t({k — /£_) luK^ + (/€+ — k) In/t-)) (12) 

However, Eq.(|T2]) ignores the [0,m] vertices: Although Eq.(|^) is valid in this case also, 

m + 1 X 

Xq = 1 X, Xm = — ; m > K, 

m m 

the energy expression of Eq.(^) differs from the Eq.(|ll]) form: 

/{0,m} = — (1 — x)(u — rlnm); m>K (13) 

Formally, /{0,m} as a function of a real variable {m-^ fi) would achieve its minimum at 
fi = exp(u/r + 1), hence, one of the nearest integers to /x, /i_ and yU+, is the true candidate 
to be the minimum point. To analyse what kind of oxygen arrangement in CFs is favorable 
we should compare two competing minima of Eqs.(p!2D-(p!3D. Such an analysis is put into 
Appendix A and it says that the global minimum of the free energy / is situated on the 



[0, fi±\ vertex if k+ < /i and on the [k_, «+] vertex if «;_ > /i. Switching over the vertices 
occurs at some fi between /€_ and k+. 

Less rich is the information about "vacancy fragments" (VFs), ie, those pieces of the 
chains which are free of oxygen atoms. Nevertheless, one can define the average length of 
such VFs and express it in terms of the average length of CFs, J, as 



fc=i 

Note, that there is a broad distribution of VFs over their lengths around /i/k. 

The conclusions of this subsection are applicable to the ensemble of non-interacting chain 
(and vacancy) fragments: 

• If the oxygen content, x, is not close to 1, ie, in the moderate-zt case, the preferable 
length of CFs is determined by /i. The average length of VFs could be estimated as 

mi,~/i//t. 

• If /«+ exceeds fi, the average length of CFs is determined by n. They are separated 
by single oxygen vacancies. 

In the following sections the role of the configurational entropy will be revealed and the first 
conclusion of this subsection will be significantly revised. 

B. The role of configurational entropy. 

By the next step we incorporate the configurational entropy term into the free energy of 
a Id chain. This will be done in order to estimate the average length of chain fragments, 
which is available in the neutron and X-ray diffraction studies (see |l]-§[). A theoretical 
prerequisite is collected in Appendix B where the "probability" set, {xm}, is derived (recall, 
that the formulas below stand upon Eq.(|]) (see also p^ , |25[| )): 



e + 2^ ^Jm{m + l)z 



^ 

e^ 



^0 = .^ ^ ^ ^. ^,. m (16) 

e + 2_^ y/m{m + Ijz 



m=l 



Certainly, constraint (||) is satisfied by Eqs.(|l^)-(0) where $ = 02/^- z plays the role of 
fugacity and could be determined from Eq.(^, rewritten in the following form: 

8 



e* + 51 v^^" 



g = Ei^l (Ij) 

We perform the summation over m by using substitution z = exp—l/2mQ, assuming that 
mo^l. Then, for example, in the main order one obtains 

roo 

J2 Vmrn^'z"' -^ y^^r^ / dxy/xx^'e"'/^ = V2^i^m]^+\2n + 1)!! (18) 

m=l ■'^ 

One may check by using the Poisson summation formula that the main order correction to 
Eq.(p!8D measured in the Eq.(|18D units is expressed through Riemann's zeta function, 

2 fn 3 \ / 3 

■ cos — TT + tTT C h^ + o 



(47rmo)"+3/2 V2 4 / ^ V 2 

and small at any rriQ of the real interest. 

Following such a summation recipe we get from Eg . ([l7|) : 

e* + v^27rmn 



e* + 3V2^ml^^ + y^m-^J^ 



Eq.(|l9l) can be properly solved by substitution: 

-S 

S 3V271 



5/2 1 - -S" e* 
^0 = o o /TT- (20) 



which assumes ttiq to be neglected as compared to e* in the r.h.s. of Eq . ([T9|) . On the 
contrary, when 5'<^1, e*<^mg' and the analogue of Eq. f^Of ) reads 



For the average length of CFs, i, one assumes the following definition: 



k=l 1 — S 



Y,Nk S -xo 
fc=i 



(22) 



One may easily check that 



3mo, mo > 1 (23) 



independently of the 5'-range. The same estimate can be found after minimization of mXm 
over m. However, such a distribution is very broad, its half-width exceeds 2i. This circum- 
stance makes the situation significantly distinct of that discussed in subsection II. A. 
The temperature dependence of i is exponential when S is not very close to 1: 

9 



. _ C. \ 2/5 

-\ e'"/' (24) 



and practically T-independent 

1 = 1/S (25) 

for S" in a close vicinity of 1. 

By using Eqs. (|I^ , (|TB|) and (^) the average length of VFs can be estimated as 

4 = j^ = 1 + -=^ ^ 1 + 3mo-^ (26) 

^ — ^0 V 2vrmo i — o 

At 5*^1 £t,~l, otherwise £v~IS/{l — S). This conclusion, concerning VFs, is similar to 
that of subsection II. A. 

C. In-chain structure factor and density-density correlation functions. 

Here we derive the equations and suggest the computational recipes of the in-chain 
correlations. Apart from the other values they allow to compute the in-chain irreducible 
averages which enter the high-temperature expansion of the free energy. 

We start with a calculation of the Fourier components of the oxygen-oxygen correlation 
function, expressed through the distribution of CFs over their lengths in the following form: 

B{q) = sl((^l-^Y.^ x^e^''^-^^) j ' + cc.^l - 1 1 (27) 

In order to prove Eq.([77|) one reminds a convention of the chain fragment definition: If 
{mi,m2,m^, . . .} is the sequence of lengths of the CFs, then the coordinates of oxygen 
vacant sites are {mi + 1, mi + m2 + 2, rrii + m,2 + m,^ + 3, . . .}. The Fourier component of 
the oxygen density (not averaged!) is simply 

g(q) = J2 e^"" - {e*'?('"^+^) -|- e^«('"i+'"2+2) + . . .} (28) 

n=l 

The first sum in the r.h.s. of Eq.(p8|) is zero for any non-zero q and we omit it because 
the subject of interest is the irreducible averages. The probability to find a CF of the 
length m in the chain fragment ensemble is Xm/ J2k=o^k = Xm/S. Averaging g{q) over the 
distribution of the CF lengths according to the auxilary formula 



1 

S 



{e'n = e E e'^'^a:. 



m=0 

one obtains: 
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{g{q)) = - ^|e*9(™i+l) ^ ^^gi„^,+m2+2) + . . . j^ 



'-' rn=0 



iq{m+l) 






1 

's 



m=0 



Jq{m+1) 



Xr, 



According to a definition the in-chain structure factor reads 



^(^) = 77 ( E ^n 



M 



I tfi^ O 



i(3(T-i-r2) 



\?'l,»-2 



M 



(29) 



^{9i-q)9iq)) = ({e^«(™^+^) + e^''(™^+™^+2) + . . .} 



To evaluate Eq. (|30|) we introduce tlie auxilary function: 



(ff ^ / ( J.q{mi+1) , iq{mi+m2+2) , , iq{mi+m2 + ...+mk+k)\ 



^ r -i<7(mi+l) _|_ -iq{mi+m2+2) , , -iq{mi+m2+...+mk+k)^\ 



(30) 



(31) 



Let us define ^ = (^e*^'-™"^^-'). Evidently, |^| < 1 for not a regular distribution of CFs in a 
chain. Now $fc can be rearranged as 

$fc = (i + r + • • • + ic)'-') + (e + 1 + r + • • • + (D'"') + • • • 



fc— 1 I /(^*\A;— 1^ 



^/c-l^^/o-2^___^^^^^ ^^(^_^)(^^^*)^^_^^(^/c-l^(^ 



'l-^k 



+ 



-e + A;e' -{k- i)e+^ 



+ c.c. 



yl-e 2J ■ (1-0^ 

Going to the fc — > oo limit according to the rule lim,^'^ ^oo k/Afe = S we obtain 



B{q) = S 



1 



1 



1 - e 1 - e 



1 



(32) 



which is simultaneously Eq.(^). 

Value B{q = 0) is a matter of interest, because it enters the thermally induced interaction 
terms of the high-temperature expansion of the second order (see Appendix B). A small-g 
expansion of B{q) yields a direct reminiscence to the Lorenzian shape of the structure factor: 



B{q) = S 



^X(2) 



l + f^(3(x(^))^-4.(3) + ^) 



(33) 
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where 

m=0 

Hence, 

B{q ^ 0) ^ S{Sx^^^ - 1). (34) 

This expression can be simply evaluated at 5*^1. Actually, according to Eqs. (p^Sf) , (|20|, (pT]) 

(2) _ 5!!v^2™^ _ _ 5 

In this case B{q ^ 0) ^ 2S/3. It tends to increase with increasing S: 

B{q ^ 0) -^ 5\(2) ^ 5^2(1 - 5)mo (35) 

where mo is supposed to be large and temperature dependent (see Eq. (|20|) ). 

The prefactor in the g^-term of Eq.(^) can be interpreted as the square of the in-chain 
correlation length, A^, (positive or negative!). So, after doing a simple, but tedious algebraic 
exercise, one obtains: 



lA^ 



f(l-5)^-f(l-5) + | 



ml (36) 



The interesting feature of Eqs. (|33D and (|36D is that the central peak at g = sphts when 
l-S> (14- v^)/15 ^0.757. 

Evidently, the knowledge of B{q) allows to perform estimation of the irreducible density- 
density correlation function: 

{{nraUra+r)) = {n^rim+r) " {n^Y = / 7rB{q)e"^'' (37) 

Note that the important sum, a{q) = S^^e^'^ J2m=o ^m^'''^"^ , entering B{q), can be asymptot- 
ically evaluated in the case of oxygen- rich chains (ttiq ^1). It can be transformed in the 
Eq.(p!8D manner: 



.,lW5e-™|!(l^||™,r!^_ (38) 

1 + V2^e-'^mo 



so, the general expression of B{q) takes the form: 

^/ X ^ g^ 2(l-Mcos^) + a(l-M^) 

2(1 — cosg)(l-|-Q;(H-MCOS'0)) + 2Ma sin ■?/' sin g + a'^{l +u'^ —2ucos{q+^{j)) 



where a = \phxmQ e~*, m = (1 -|- 4g^mo)~^^^ and ^ = |arctan2gmo- One notes that by 
using Eq.(pOD all the parameters, entering Eq.(^), can be expressed through mo- Namely, 
a = (1 — S)/{SmQ) and small. In the 5*— >-0-limit a is no more small, on the contrary, a^ 1 
and B{q) takes the form: 

1 -M^ 
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III. STRUCTURE FACTOR OF INTERACTING CFs. 

In the preceding section we were interested in the behavior of the chain fragment en- 
semble: First, in a hypothetical situation provided a thermal equilibrium of oxygen in a 
basal plane is available at rather low temperatures, second, in more realistic case when the 
configurational entropy is not neglected. The results obtained in the latter case seem to be 
reasonable as regards the in-chain correlation functions, although the interaction between 
the chains has been ignored. However, the role of such an interaction is important in many 
physical quantities which are sensitive to distinctions between Ortho-I and Ortho-II phases, 
such as concentration of holes leaving CFs, nh, concentration of paramagnetically active 
CFs, ripm, and so on. In this section we derive the formulas, describing that part of diffuse 
scatterring caused directly by oxygen atoms and vacancies. We apply these formulas in 
numerical calculations, and visualize them in Section IV. 

Treating the problem in terms of the set {xk}, we outline in Appendix B the method 
which incorporates the interaction between nearest chains into some scheme of statistical 
mechanics: This allows to determine the correlated set {x^}, depending of the type of 
orthorhombic phase. The goal of this section is to develop a similar approach for correlation 
functions, in particular, for the irreducible part of the oxygen density-density correlation 
function. The latter is associated with the intensities of diffuse scatterring due to a partial 
disorder in the CF ensemble: 

^(q) = ^EG'(ri,r2)e^^(r-r^) (41) 

-'*' ri,r2 



where 



G'(ri,r2) = {{rir^rir,)) = {rir^rir^) - (^ri)(nr2) (42) 



and A/" = A/'^A/'ch- It is evident, that for ri and r2, both belonging to the same chain, the 
main order contribution into G{ri, r2) relates to the in-chain correlation function discussed in 
subsection II. C. Provided ri and r2 are separated by n chains, the main order contribution 
appears to be 

G(a,/xo;« + n,/x„)= --^ ^... E (Ko<o)) («'<f )) " " " ((O^D) (43) 

^ -^ ^ mo run-i 



or. Eg . (|43|) can be read as 

Go.A%) = jf EG'(a,/io;a + ri,/in)e*'^''('^«-^")=f-^j 5«(g,)5,+i(g,) . . . 5„+„(g,) (44) 



As regards the Eq.(|i^ derivation we refer the reader to Appendix C. We use in Eq.(PD the 
y component of the wave vector following the convention about predominant orientation of 
chains along y-axis. 
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Let us perform summation in B{q): For any separation n between the chains the leading 
term of the order (V3/T)" will be only kept. 

The Ortho-I case is rather simple: All the elements B^ = B(qy) in Eq.(^) are a- 
independent. Summation over all the elements Gn{qy) in Eq.(^TD yields: 

S(q) = i? (1 + 2 ^ cosng. (-^^)"| = ^ V + XT ^ \ (^5) 

I n=o ^ ^ ^ J l^i + ^e^''-^ 1 + — e-^^-5 

In the case of the Ortho-II phase the analogous expression is more complex {Ba alters from 
chain to chain, taking the values, say, Bi{qy) and B2{qy) for oxygen-rich and oxygen-poor 
chains, respectively): 

^(q) = ] ^^^ - 2^^i52 cosg, + + (^^) (Bi + B2)BiB2 cos2g. 



-2^^) BfBlcos3q, + ... 



^ T — + C.C. 






1 - -^e'^'^-'BiB2 



Bi + B: 



2 



(46) 



Anticipating the results of numerical calculations let us find the B{q) shape in the Ortho- 
I phase in the vicinity of the Ortho-II Bragg peak (^ = Sq^, Sqy): This point of the reciprocal 
space is a suitable candidate to be used in determination of correlation lengths along and 



across the chains. One simply checks that Eq. (|45|) can be written as 



^ _ fVsB{6qy 



B{q) = Bi5qy) ^ '- -^ (47) 

1-2 jr^cos6q^+ I jr^j 

The correlation length across the chains is associated with the denominator in Eq.(^): 



' (1 



VsBjO) 
T 



SO, A^ is positively defined. The A^ expression can be also easily obtained from Eq . (^71) : 



VsBjO) 



For A^ see Eq.(|6|). 

If V3-B(0)/T^l, ie, provided the high-temperature expansion holds, A^. is of the order 
of the lattice constant, whereas A^, ~ |A| ~ uiq. 
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IV. NUMERICAL CALCULATIONS. 

Below some analytical results of the preceding section are visualized after numerical 
calculations. In order to make a selection of the fitting parameters for the theory presented 
above we refer to the plateau-like behavior of rih vs x which should imitate the 60 K plateau 
of Tc vs X in YBCO. In Fig. 3 the n^ix) dependence is shown at V2 = V3 = 1 and Tg, varying 
with interval 0.08 from 0.46 to 0.86, which seems a reasonable choice for interpretation 
of experimental facts. The regions of stabilities of Ortho-I and Ortho-II phases can be 
clearly distinguished in the set of curves of Fig.3. The curves, illustrating the monovalent 
copper amount vs x are shown in Fig.4. A propos, monovalent copper atoms occupy 2- 
fold coordinated sites, which are out of CFs. This quantity, 'n.2_foid, is certainly not in a 
contradiction with the experimental data of Refs pT] , P^ . Closely related to n2-foid are the 
concentrations of copper ions in 3- and 4-fold coordinated positions, ie, copper atoms at the 



ends of CFs and within them, respectively (see Figs. 5, 6). According to |^ the information 
about distribution of copper atoms in differently coordinated positions can be obtained 
in NQR experiments. The quantity, which also reflects the structural peculiarities of the 
Ortho-I and Ortho-II phases, is the average length of chain fragments shown in Fig. 7. The 
last physical quantity of the reference set available for checking experimentally could be the 
concentration of paramagnetically active chain fragments, npm. The theoretical curves have 
typical forms shown in Fig.8. The experimental research program, employing the polarized 
neutron technique is in the beginning (see Ref. [^). Also the concentration of paramagnetic 
CFs at low temperatures could be measured due to the Shottky effect. Below we give shortly 
the explanations of how the afore-said concentrations can be evaluated in the frame of the 
generalized lattice-gas model. 

A. Concentration of holes quitting the chains. 

The general expression of the concentration of holes transferred from chains has a form: 

1 " 

2 n=l m=0 

Factor | enters the equation because one oxygen deficient plane supplies holes into two Cu02 
planes. Function Pr(?T,, m) is the probability for a chain fragment of the length n to have m 
0^~ ions within and can be obtained in the frame of some strongly correlated model. Here 
we accept a simplified consideration according to which the most probable configuration has 
to be substituted Eq. (|50[) . For very long chains the most probable m at fixed n is simply 
^n, where ^ = lim„^oo^n- Although for shorter chain fragments the relationship between rn 
and n is more delicate, it could be approximately expressed as 
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m{n) ~ nmt(^n) (51) 

where function nint denotes the nearest integer number. So, in this case we rewrite Eq.(pOD 
as follows: 

B. Concentration of paramagnetically active CFs. 

Again, starting with some microscopical model for a finite CF characterized by n and 
771, one obtains, in principle, the average value of S'^ . To simplify the scheme (see also Ref. 
P4|| ) we first recall that there is no charge transfer from short chain clusters with one or 
two oxygen atoms, both are spin-singlet fragments. For a cluster of three oxygen atoms, 
n = 3, a spin-1/2 state forms, contributing to the paramagnetic Curie constant. On the 
other hand, for n = 4, one expects again a singlet state, because the most probable is the 
m = 3 state. In the longer clusters, n > 5, singlet and doublet states occur with roughly the 
same probability, at sufficiently high temperatures. So, the paramagnetic Curie constant 
is proportional to the concentration of paramagnetically active CFs, ripm, which may be 
approximated by 

n^^ = {xz)^ -^{Xn) (52) 

^ n=5 



C. Cu+ amount. 



The experiments pi| , |22| were very important for understanding of the charge transfer 
mechanism in YBCO. Because only those copper site which are surrounded by oxygen va- 
cancies give a contribution into the monovalent copper amount, their concentration coincides 
with (xo). 

D. Average lengths of CFs and concentrations of Cu atoms in 3- and 4-fold 

coordinations. 

Acceptable in applications to the Ortho-II phase are the formulas given below and ex- 
pressed by means of the Xq values: 

Ac) 

4 = — TT (53) 

where (c) relates to either rich (r) or poor (p) chains; their average characteristic takes a 
form: 
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Noting that 2(A^i + A^2 + • • •) is the amount of copper in the 3-fold coordination in a chain, 
one easily obtains n3_foid 

ns-foid = 2 -2x-2(a;o)- (55) 

^4-fold can be simply derived due to a sum rule: 

n4-foid = 1 - ^2-fold - ^3-fold = 2a; - 1 + {xo). (56) 

E. Structure factor: Diffuse scatterring on the CF ensemble. 

Figs.9(a-f) illustrate a diffuse part of the structure factor calculated numerically accord- 
ing to the scheme developed in this paper. Remind, that this part is due to structural 
imperfections in the CuOa: plane. The fitting parameters, V2 and V3, are common with the 
values used for the plots in Figs. 3-8. Calculations were performed at Tq = 0.62: Figs. 3-8 
also include the curves relating to this temperature. 

The following features of the structure factor shape are necessary to be mentioned: 

• Its value enhances on the Brillouin zone boundary, (±1/2, g^^) (in units of 2tt), it 
reflects a tendency for a system to be arranged as a double period structure on a short 
wavelength scale. This fact would become especially pronounced if the Ortho-I phase 
were forced to substitute a true Ortho-II phase - in this case a short-range oxygen 
atom - oxygen vacancy correlation across the chain should be significant. 

• The Ortho-I phase exhibits a clear down-fall behavior of the structure factor in the 
(±1/2, 0) points of the reciprocal space. This reminds a central peak splitting discussed 
in subsection II. C. 

• The maximum diffuse scattering intensity is displayed in the vicinity of the Ortho-I- 
Ortho-II phase transition. 

V. CONCLUSIVE REMARKS AND SUMMARY. 

This paper is devoted to the correlations in the oxygen deficient planes which are the im- 
portant structural elements of HT superconductors, YBa2Cu306+2; and related compounds. 

There were many experiments performed to reveal the fascinating properties of such 
compounds. Some of them are quoted in this paper in connection with superstructure 
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formations, charge transfer mechanism, magnetic behavior of CFs. There were also many 
attempts to describe CuO^ planes theoretically. Many of them invoked into consideration 
the oxygen subsystem only. We already mentioned the conventional lattice-gas approach 
PJ6|JT5Hr8| ; the approach based on the screened Coulomb potential |Q has been also applied 
to the oxygen subsystem in order to explain some hypothetical superstructures. The model 
utilized in this paper looks as generalization of the lattice-gas model (see also |]2^), but the 
role of copper in formation of alternating Cu-0 CFs and charge transfer mechanism is not 
ignored. 

The statement, utilized in this paper, concerns the copper-oxygen chain fragments: Many 
phenomena in YBCO are caused by them. The function of the CF ensemble to control the 
YBCO behavior is closely connected with the charge transfer which could be conditionally 
subdivided into two parts: the intrachain charge redistribution and that one from chains 
onto Cu02 planes. The former is the oxydation of monovalent Cu ions, forming the CFs, to 
the bivalent state According to this process almost all oxygen atoms in not very short Cu-0 
CFs would be subjected to transformation into the monovalent state where it not for the 
hole transfer to outside: Around 30% of oxygen atoms in chains persist their monovalency, 
the rest oxygen constituents appear to be in the bivalent state [|l9l. It is noteworthy, that 



the distribution of oxygen valencies is neither frozen, nor ordered in any charge density 
superstructure, but due to strong p-d hybridization oxygen holes form strongly correlated 
systems of finite lengths as discussed in |^ . 



According to the detailed investigation by Jorgensen et al |^3| of the orthorhombically 
arranged CuO^: planes the concentration of oxygen atoms in "wrong" 0(5) positions is 
negligibly small up to rather high temperature ~ 600°C. The situation reminds the soft 
matter: Cu-0 fragments behave as long molecules strongly oriented and weakly correlated 
to each other in chains. 

Any oxygen migration stops approximately at and below room temperature. So, this 
circumstance does not allow to appear those CF structures whose existence could be a 
result of a simple energy competition. Actually, distribution described by Eq.(|15D has a 
maximum at itiq, but its width is of the order 2mo, hence, the role of the entropy becomes 
decisive. 

Below we resume the method and results presented in the paper. 

• The principal variables over which the statistical averaging can be done are the prob- 
abilities, {x^}. 

• The broad distribution of CFs around Sttiq (see Eq.(^)) can be utilized for calculation 
of such quantities as the hole transfer from chains, concentration of paramagnetically 
active chain fragments, the monovalent copper amount, etc. They evidently depend 
on the oxygen content, x, and quenching tempearature, Tq. 



• The in-chain structure factor can be used for interpretation of the diffuse scatterring 
results in the Ortho-I phase. Sphtting of the Id central peak could be a counterpart 
of the analogous splitting of the diffuse maximum on the Brillouin cell boundary, 

(±1/2,0). 

• Oxygen-oxygen correlations in the basal plane are illustrated in Figs 9. The maximum 
of diffuse scatterring corresponds to the Ortho-II-to-Ortho-I transition concentration. 

• Measurements of the correlation length along the chains can be employed for experi- 
mental determination of effective coupling constants. 

The message of this work is clearly addressed to those experimentalists whose field of interests 
lies in atomic and magnetic correlations in YBCO and related compounds. Neutron and 
X-ray studies seem to be the most appropriate technique for this. Nevertheless, the method 
developed in the paper is possible for evaluation of thermodynamical quantities, such as, 
isothermical susceptibility, Shottky anomaly, etc. 
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Appendix A. 

Taking tlie formal derivatives of f{i,j} (see Eq.(|TT|)) over i and j we get 



dj (j - i] 



^^^'''^- -r(l-x)/:i^p-^-ln^-) (A.2) 



di (j — 2)2 

Substituting z = i/j < 1 we notice that the sign of derivative ( |A.1| ) coincides with the sign of 
expression z — l — lnz and positive. Let i be fixed, hence, f{i,j} approaches the minimum 
as a function of j at the greater nearest integer to k, say, k^. The same substitution makes 
the sign of derivative ( [A .21 ) identical to the sign of expression 1 — lnz — l/z which is negative. 
The minimum of f{i,j} as function of i is situated at the lower nearest integer to k, say, 
/€_. Hence, f{i,j} attains the global minimum on the manifold of vertices, defined by 
Eq.(^ with < z < j, on unique vertex [/€_, k^]. 

In order to search, what kind of vertex, [/?_ , k_|_] or [0,m], will be realized as a proper 
one, we start with the case /€_ >/i = exp(l + u/r). Functions 

fi = u — t((k — K)lnK^ + (k,^ — K)lnK,) and /2 = — (u — Tlnm) 

m 



differ from those of Eqs.(|l^) and (|T3D by a common multiplayer. Because in this case m = k+ 
is optimal for /2 (see Eq.(|TB|), 

/i — /2 = M 1 — T \ {k — kJ) In «;+ + (/?+ — k) In K- In /«+ 

V ^^+J \ /«+ / 

= r (In/i — 1 — /t_|_lnK_ — /t+ lnK4.) (A. 3) 



appears to be negative and favors the [«;_, /?+] vertex. When deriving Eq. (|A.3|) we used the 
relation In/x = u/t + 1. Contraversially, if k+ < /x we substitute /2 by /2: 



Then 



/2 -^ /2 = — (u - r In K+) > /2 



/1-/2 >fi-f2 = r— -{\nfi~ 1 - K+\nK_ - K+lnK+) >0 (A.4) 



To complete the proof we consider the last possibihty, k_ <yU< k_|_. Again k+ substitutes 
for fi and the sign of (/i — /2) coincides with the r.h.s. of Eq. ([A.3|) which turns to zero at 



1 5 

^ 2 24k_ ' 
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Appendix B. 

Deriving the equations of statistical mechanics for oxygen in a deficient plane we suppose: 

• A negligible probability for oxygen to occupy a position across the chains; 

• The largest energy contribution comes from the interchain degrees of freedom; 

• Apart from the entropy contribution to the total energy of the CF ensemble we take 
into account the oxygen-oxygen repulsion on the nearest chains, like in the lattice-gas 
model (cf, m&- 

• The constraint used is the total oxygen content, x (see Eq.(^). 

As regards the second supposion it denotes a preferable role of strong fermionic correlations 
within . . .-Cu-0-. . . chain fragments as compared with interchain interactions which are 
effectively screened in the metallic state of YBCO. The oxygen ordering in deficient planes 
reminds the soft matter case, e.g., the nematic liquid crystal. Whereas the CFs weakly 
correlate to each other they are strongly coupled from inside and preferably oriented in 
chains. 

Let {Nl^, A", N2, . . . , N^, . . .} be the total numbers of CFs of various lengths (for their 
definitions see the beginning of Section II) in the a-th chain. A complete number of different 
configurations of all those CFs has a form usual for a mixture of ideal gases: 

N(}\N^\N^\ . . . N^\ . . . ^ ■ ^ 

Before giving a definition of the partition function of the CF ensemble in the oxygen deficient 
plane let us omit for a while the Va-term and consider the auxilary partition function, Zq, of 
the a-th chain with two constraints imposed: (|^) and (^. Playing the role of the Lagrange 
multiplayers are the chemical potentials, fi"" and /ig, conjugated to the above mentioned 
constraints. With such a prerequisite Zq takes a form: 



Zq = J2 expM < ^a In ^a - Z! ^m Inx^ - 7^ X! 

{x"} I rn=0 -^ m=l 



m(j){rn)x'^ 



/^ V-.^ , lN.a /^O 



+ ^ E (^ + 1)^::. -h^S-\ = Y: exp A-^dx^) (B.2) 

^ m=0 ^ J {xD 

where summation runs over the possible distribution of CF lengths, y4"({x^}) is the corre- 
sponding "action". 
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The equilibrium distribution of "probabilities" {x^} achieves at the niaxinial "action" 
and can be expressed as follows: 

m=0 



E 

m=0 



e 



-m{4>{Tn)-fi°')/T 



qa _ rn=u /-p a\ 



^ (m + l)e 

m=0 



-m{(p{m)-ti°')/T 



In principle, /Xq can be also expressed through /x": 



m=0 



So, any Zq can be considered as a function of /Iq,. Substituting (f){k) in Eqs.( |B.3|) -( p.4| ) by 
its approximation of Eq.(|^) one obtains functions ([T^)-(]T7D. 
Let us include the Vs-term into consideration: 

Z= Yl e^"(K})(-."+i; (°."+i = exp I -^ Y. «^M (B.6) 

{x-} V ^ m+1 J 

Then it is suitable to expand (^°'"+^ over the cumulants. Such a scheme is well-known and 
we only outline it below {r] = Ha C°'"+^): 

K} K) V ^ / V ^ / 

where < . . . >o an average over the ensemble of non-interacting chains, Zq = Ha 2^o- For 
example. 



-^ a 

The expectation values, entering the (?7)o-term of the high-temperature expansion, can be 
expressed through the S'-values, supposed to be known: 



^ n-n-+i = ^7^ n°+i = 5^(1 - 5")(1 - ^"+^) = (2s - l)A4h + J] ^"^" 



+1 



(?7^)o consists of the expectation values of the forth order, nfn1~^ ''^j'^j i "which completely 
decouple if/3>a + l oi (3 < a — 1. Also, 

<-^<n^"nJ+^ ^ ^^^^ ^^^+T (7F^2 ^ (f^nyj)) ) (B.7) 
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where 

(«<))o = <<-^^ (B.8) 



According to a definition ((^f^j))o = at a 7^ /3. Note, that Eq.( |.B.8D is invariant under 
the oxygen atom <^ oxygen vacancy transformation (n^ ^^ 1 — n"J. 
There is another non-trivial contribution of the r/^-term: 

«+^n^"nJ+i ^ (n^^ + («np)o) (^^^^' + («+'nJ+i))o) (B.9) 

All the decoupled terms, coming from all the orders of the high-temperature expansion, can 
be recollected into 

exp(-^E^"^"^')- (B.IO) 

The terms ((nfn"))on"~^ 72"+^ (see Eq.( [B.7| )) give rise a thermally induced interaction of 
the next-to-nearest chains. The recollected contribution of Eq. (p.lO| ) is of the main order in 
the V3/T expansion: It formally corresponds to the mean- field approach. The Vs-exponential 
term (Eq. (|B.lCI| )) simply contributes Zq and does not formally influence on Eqs.( |B.3|) - ([B.4| ), 
but Eq.( |B.5| ) changes its form 

■m=0 ^ ^ 



It is noteworthy, that Eq.( |B.ll"D is a system of non-linear equations which couple the 



densities of oxygen vacancies in nearest chains. Actually, S'" = ^"(/i") (see Eq. (|B.4|) ), 
hence, the r.h.s. of Eq.( [B.ll] ) is some function of S*""^, 5*°, 5"^^, whereas the l.h.s. of 
Eqs.( |B.lTD does not depend of index a and contains the function of the global chemical 
potential, /iq, conjugated to the (^/ofta/ density of oxygen atoms (vacancies). 

Appendix C. 

The correlation function entering the r.h.s. of Eq.(^2p can be expressed through the partition 
function of non-interacting chains, Zq, as follows: 

Y, ZoUr^rir, exp -^ E«+^ 



E^oexpl-^5:«+^ 
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V.^. „, .,^,, ^2 



{nr,nr,)o - -^ 5I(^ri?^r2«"*"^)o + T^;^ im(^ri^r2«+^^?^f ^^)o + 






In order to deal with irreducible averages one needs the (rir) definition: 



(C.l) 



(nr) = '-^ j^ "-^ 

{x-fe} \ a,i 



V.^, ^. „,^,. K 



2 



)o - f EK<nr^)o + ^ E E(^r<nr^<nf ^)o + . . . 

1 - f E(«^')o + ^ E E(«^^«"'')o + • • • 



(C.2) 



It has been known since the beginning of sixties [^ that in Ising-like systems a corre- 
lation function can be expressed as a power series over irreducible averages, the so-called 
cumulants. A true expansion must only contain connected diagrams. In our particular 
problem the interaction involves atoms of nearest chains, hence, a series expansion of the 
irreducible average, ((n^n"'*'"')), starts with a term proportional to (Vs/T)"': 

{{n>rn) = 7^) E((«))oE(«t'<'))o---E((Or"))o (C.3) 

^ -^ ' mi 7712 nin 
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Figure Captions. 

Figure 1. The structural phase diagram is shown qualitatively. The bottom is confined 
by room temperature, Tr: Below Tr oxygen is practically frozen. Probably the Tetra-to- 
Ortho-II phase transition is of the 1st order and the T — x phase diagram contains the phase 
separation area. Hypothetically, the superstructures of longer periodicities, like Ortho-Ill, 
could appear between Ortho-II and Ortho-II areas. 

Figure 2. The lattice-gas model definition, involving the oxygen-oxygen interactions only 
(circles and crosses are for oxygen and copper atoms, respectively). 

Figure 3. Uh vs x: Six curves are shown for the equidistant set of quenching temperatures, 
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from 0.46 to 0.86; V2 = V3 = 1; The higher temperature, the smaller nh- 

Figure 4. ncu+ (^2-foid) vs x: Six curves are shown for the equidistant set of quenching 
temperatures, from 0.46 to 0.86; V2 = V3 = 1; The higher temperature, the smaller n2-foid- 

Figure 5. n3_foid (copper in 3-fold coordination) vs x: Six curves are shown for the equidis- 
tant set of quenching temperatures, from 0.46 to 0.86; V2 = V3 = 1; The higher temperature, 
the larger n3_foid- 

Figure 6. n4_foid (copper in 4-fold coordination) vs x: Six curves are shown for the equidis- 
tant set of quenching temperatures, from 0.46 to 0.86; V2 = V3 = 1; The higher temperature, 
the smaller n4_foid- 

Figure 7. Average length of CFs: d (solid curves), dr (dashed curves), dp (dotted curves); 
V2 = V3 = 1; Tq = 0.46, 0.62, 0.78. The higher temperature, the smaller ts. 

Figure 8. ripm vs x: Six curves are shown for the equidistant set of quenching temperatures, 
from 0.46 to 0.86; V2 = V3 = 1; The higher temperature, the larger Upm- 

Figure 9. The structure factor, i3(q), as defined in the text. V2 = V3 = 1; Tg = 0.62; 
Oxygen content changes equidistantly from x = 0.45 (a) to 0.95 (f). Chain orientation is 
supposed along the y-a.xis. 
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